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A TEST FOR MONOMIAL CONTAINMENT 


SIMON KEICHER AND THOMAS KREMER 


Abstract. We present an algorithm to decide whether a given ideal in the 
polynomial ring contains a monomial without using Grobner bases, factoriza¬ 
tion or sub-resultant computations. 


1. Introduction 


Let 1 be a field. Given an ideal I C K[Ti,... , T r ], the monomial containment 
problem is to decide whether I contains a monomial. Equivalently, one is interested 
in whether the intersection V{I) flT r of the zero set V(I) C K with the algebraic 
torus T r := (IK ) r is empty. The monomial containment problem occurs frequently 
when determining tropical varieties jl] or when determining GIT-fans [10]. The 
usual approach is via Grobner bases: I contains a monomial if and only if the 
saturation I : (Tf • • • T r )°° contains 1 £ K[Ti,..., T r \. This can also be decided by 
a radical membership test: I contains a monomial if and only if T\ ■ ■ ■ T r £ \fl. 

In the present paper, we provide a direct approach involving neither Grobner 
basis computations nor (sub-)resultants or factorization of polynomials. We con¬ 
sider more generally the following problem: given a polynomial g £ 1K[7i,..., T r \. 

- 7 * 

prove or disprove the existence of an element a; £ IK such that 


(1) f(x) = 0 for all / £ I, g(x) ± 0. 


Clearly, setting g := T\ ■ ■ ■ T r £ K[Ti,...,T r ] in Q, the existence of such x is 
equivalent to the monomial containment problem. Our algorithm, Algorithm |4.1[ 
proceeds in three steps: 

(i) Compute finite subsets Si,...,S m C K[Ti,...,T r ] that are in triangu¬ 
lar shape and polynomials gi,, g m such that the solutions of (JTJ) are 
preserved, i.e., the zero sets satisfy 

v(i)\v( g ) = |jm)\n<?*) c r. 

(ii) Making certain variables Tj invertible, we obtain a function field L and 
an embedding l : K[Ti,..., T r \ —► L [T ^,..., Tk s } such that the embedded 
equations i(Si) are dense , i.e., each variable Tl. corresponds to an equa¬ 
tion. 

(iii) Then an element x £ K satisfying 0 exists if and only if the minimal 
polynomial of the class i(gi) £ L[Tfc 15 ... ,Tk s ]/(Si) is not a monomial for 
some i. 


Experiments with our implementation of Algorithm 4.1 suggest that it is compet¬ 
itive for certain classes of input; for instance, it usually beats the Grobner basis 
approach when a solution exists, i.e., the ideal is monomial-free. 

Note that the idea behind step (i) of the algorithm is quite common and similar 
concepts have been used by several authors for a more explicit study or even the 
explicit computation of solutions. See, e.g., □Q GO El 0H7| for a series of papers with 


2010 Mathematics Subject Classification. 13P05, 13P10, 13P15, 14Q99 . 
The first author was supported by the DFG Priority Program SPP 1489. 


1 





2 


S. KEICHER AND T. KREMER 


Grobner basis-free algorithms for systems of equations. The methods of Wang m, 
Thomas 031CH] as well as Bachler, Gerdt, Lange-Hegermann and Robertz 13] can 
also deal with systems of equations and inequalities. They determine the solutions 
of such systems by means of certain triangular sets called simple systems ; their 
computation involves sub-resultant computations. All algorithms, including ours in 
step (i), share the concept of triangular sets , certain finite subsets Si C K[Ti,..., T r \ 
such that V(I) = (J V(Si) holds. The Si then give insight into the structure of the 
solution set V(I) C K . As we are only interested in solvability of 0 , we will 
only need triangular sets with weaker properties but which can be computed more 
efficiently. 

The structure of this paper is as follows. In Section [2] we show how to decom¬ 
pose the given ideal into a list of triangular sets with sufficient properties for our 
solvability test; this is step (i) in the previous list. Section [3] is devoted to steps 
(ii) and (iii), i.e., we show how to reduce the problem to a dense system over a 
function field and how to determine the solvability of such a system by means of 
minimal polynomial computations. Explicit algorithms are given in each section. 
In Section [4] we present our algorithm for the monomial containment problem. 
We compare the experimental running time of the perl implementation |12l of the 
algorithm to the Grdbner basis approach as well as to the methods of m ns]- 

This paper builds on m- We would like to thank Jurgen Hausen for helpful 
discussions. 


2. Triangular shape 

In this section, we treat item (i) of the list on page[l] i.e., we decompose a system 
as in |l]) with an ideal I C K[Ti,..., T r \ and a polynomial g £ K[Ti,..., T r \ into 
a list of finite sets of polynomials that are in triangular shape. We show how to 
compute this decomposition by iteratively applying a set of operations that do not 
change the solvability of (JT|) . 

We first define the notion of triangular shape. In the literature, they are also 
called triangular sets mmm- 

Definition 2.1. Fix the lexicographical ordering T\ > ... > T r on K[Ti,... ,T r ). 
We call polynomials /i,..., f s £ K[Ti,... ,T r \ of triangular shape if for each fj , 
there is 1 < k(fj) < r such that 

(i) we have k(fi) < ...< k{f s ), 

(ii) fj £ ..., T r ] \ K[T k ( f .) +1 ,..., T r \ holds for each 1 < j < s. 

We denote by deg T . (/) the (Ti-)degree of a polynomial / £ K[Ti ,... ,T r \ considered 
as an element of the univariate polynomial ring K[T); j ^ i\ [Tj]. Moreover, we write 

LCfc(/i)(/i) € RckUi) '■= ^[^t(/i)+lv)^r] ■ 
for the leading coefficient of the polynomial /,; considered in the ring R<k(f i )[Tk(f i )]- 

We now introduce the concept of (semi-) triangular systems. Assume / is gen¬ 
erated by polynomials j\..... f s £ K[Ti, ... ,T r ). We sort them into two sets (and 
keep track of the inequality g): polynomials that are already in triangular shape 
.F<j and remaining polynomials J-\j. 

Definition 2.2. A semi-triangular system (of equations) is a tuple (J r \j,J 7 <s ,k,g) 
consisting of finite subsets C K[T 1; ..., T r ], an integer 0 < k < r and a 

polynomial g £ K[Xi,..., T r ] such that 

(i) is of triangular shape, 

(ii) we have LC fc(/ )(/) | g for all f £l <,, 

(iii) the set {!,..., A;} contains {£;(/); / £ ^<i}, 
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(iv) for all / £ Fa and each 1 < i < k we have deg T . (/) = 0 . 

Moreover, we call a semi-triangular system (.Fb, F<^, k, g) a triangular system if 
Jb C K holds. 

Example 2.3. Define in K[Ti,T 2 ,T 3 ] the subsets .Fb := 0 and F<\ := {/i,/ 2 ,/ 3 } 
where the /* and k(fi) are 

fi := T? - (T 2 + T 3 )T u k(f r) = 1, 

h ■= t|-T 3 , fc(/ 2 ) = 2, 

h ■= T 3 2 -T 3 , *(/ 3 ) = 3. 

Then F < is of triangular shape and (Fb ,Fb,3, T) T 2 T 3 ) is a triangular system. 

Definition 2.4. A list S of semi-triangular systems is called a triangle mush. Two 

triangle mushes S and S' are equivalent if we have V(S) = V(S') with the solutions 

V(S) := U F(JbUj’ < )\y( 5 ) C It'. 

(F\jF<,k,g)eS 

For the case of a single element S = {S'}, we will use the same notions for S instead 
of S. 

Example 2.5. Consider the triangle mush S := {(Jbj, 0, 0, g)} in Kp},..., T 4 ] 
where g := T 1 T 2 T 3 and Fz 1 consists of the two polynomials 

fi := (T 3 -Ti)(T 3 -T 2 )T 2 , f 2 := (F + T 2 - T 3 )T 4 . 

—4 

Going through the different cases, one directly verifies that E(S) C K consists of 
all points {x\, x 2 , Xi, 0) and (x±, X 2 , £ 2 ,0) £ K where £ K . We will continue 
this example in |4.3[ 

Given a triangle mush S, we are interested in operations that transform S into 
an equivalent triangle mush S' that consists of triangular systems. 

Construction 2.6 (Solution-preserving operations). Let S := {(-Fn, F<\, k, g)} 
consist of a semi-triangular system. Each of the following operations produces an 
equivalent triangle mush S'. 

(i) Case-by-case analysis: If / £ K[Tfc +4 ,..., T r \ and h £ K[T 4 ,..., T r \ are 
such that g \ h and h \ fg, then one may choose 

5' := {(F n U{f},F <t k,g), (F^F^^h)} . 

(ii) Polynomial division: Consider f,h£ Fa and b £ K[T 4 ,..., T r \ with b \ g. 
Assume that for some j £ Z>o we have 

V f = ah + u, a,u £ Kpit+i,..., T r ] 

where b := LC T k+1 (h) and deg Tfc+i (u) < deg Tk+1 (h). Then we choose the 
triangle mush 

5' := {(F a \{f}U{u},F < ,k,g)}. 

(iii) Unused variable: If k < r and cleg Tfc+i (/) = 0 holds for each / £ Fu, then 
we may choose 

S' := {(.Fb,.F<i,fc-I- 1,3)}- 

(iv) Sort polynomial: If k < r holds and there is exactly one polynomial / £ Fa 
with deg Tfe+1 (/) ^ 0 and LCfc(/)(/) | g , then we may choose 

5' := {(F a \{f},F < U{f},k + l,g)}. 
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(v) Last polynomial: Assume k < r and there is exactly one polynomial / £ 
Jb with deg Tfc+i (/) 7 ^ 0. For — 1 < j < d , we write 

d j 

f = ^ ' a iTk+ 1 ; fj := ^ , 

i=0 i=0 

•?< == {/,}, -7^ := (J 7 n \{/})U{a j+ i,...,a d }. 

Then we may choose 

S' := {(J r ^,J 7 ^,k + l,gai),...,(J r (j,J r !} ] ,k + l,ga d ), 

{T □ 1 , .F <\,k + 1, g)} ■ 

Proof. One directly checks that in all cases S' is a triangle mush. For (i), each 
x £ V(<S) either satisfies f(x) = 0 and g{x) / 0 or we have /( x) ^ 0 and h \ fg 
implies h(x) ^ 0, i.e., x £ V(S'). The inclusion V(S') C V(S) is clear from g \ h. 
We come to (ii). Each x £ V(«S) satisfies 

u(x) = b(x) j f(x) — a{x)h{x ) = 0. 

For the reverse inclusion, we use b | g to obtain b(x) ^ 0. Consequently, we may 
infer f(x) = 0 from 

b{x) j f(x) = (V f){x) = a(x)h(x) + u(x) = 0 . 

Operations (iii) and (iv) are clear. For (v), we define the following triangle mushes 
for 0 < l < d: 

Si ■= {(■?□ 0 {/;}, F<,k,g)}, Vi := {(J r ^,J rJ < ,k + l,ga j y, l<j<d)}. 

Observe that by an application of operation (i), we obtain an equality of solutions 

V(S l ) = {fi},F<,k,g), (P n U {fi},F<,k,gai)}). 

As the ideal (J 7 ^ 1 U {/;}) equals (J- 1 ^ 1 U {fi-i}) and by an application of operation 
(iv), we obtain 

m) = v ({(J^- 1 U {/,_!}, T^k,g), (T^T^k + hgm)}) 

= V (Si-i U \ Vi)). 

Adding the equations stored in Vi on both sides does not change the solution set, 
i.e., V (Si U T>i) is equal to V (S;_i U T>i-i). Iteratively, we obtain V (Sd U V d ) = 
V (So U V o). Moreover, because of /o = ao and operation (iii): 

E(So) = y((J' 5 1 ,^,fc, ff )) = V((Ta\^,k + l,g)). 

We conclude that V(S) equals V(S d U V d ) = V(So U V 0 ) which in turn is the same 
as the solution set V(S'). □ 

The next algorithm transforms a triangle mush into an equivalent triangle mush 
consisting only of triangular systems. Given a triangular system (J 7 ^, !F < ,k : g), the 
idea is to reduce T^+i-degrees of an element / of the unsorted polynomials T b 
by successive polynomial divisions; afterwards, we move / into the set of sorted 
polynomials 

Given a finite set of polynomials J- C K[Ti,..., T r ], its reduction is a finite subset 
red(J') C K[Ti,..., T r ] such that 

LT(/0 | LT(/ 2 ) for all /i ,/ 2 £ red(7b), 

(LT(J-)) C (red(LT(J 7 ))), (J 7 ) = (redp 7 )) 

where we denote by LT(/) or LT (M) the leading term of a polynomial / or set of 
polynomials M with respect to the ordering defined in Section [2j Computing the 
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reduction of F means successively applying the division algorithm to the elements 
of F, see, e.g., 6 ]. 

Algorithm 2.7 (MakeTriangular). Input: a triangle mush S in K[Ti,... ,T r ]. 

• While there is S := (Jb, k, g) £ S with k < r, do: 

— Replace Fa by its reduction red(J r n)- 
— If there is f £ Fa with deg Tfc+1 (/) > 0, then: 

* If there is h £ Fa \ {/} with deg Tfc+i (ft.) > 0, then: 

• Perform a polynomial division of / by h in the univariate 
polynomial ring R := K(Tfc + 2 ,... ,T r )\T k+ i] to obtain 

/ = a'h + u' £ R. 

■ Set b := LC fe+ i (h) £ K[T k+2 , ■.■ ,T r ] and j := deg Tk+1 (h)+ 
1 £ Z>o- With a := Wo! and u := Vu' £ K[Tfe + i,... ,T r \ 
we then have 

Vf = ah + u £ K[Tfc + i,...,T r ]. 

• Redehne S := (S\ {5}) U {S', 5"} where 

S' ~ ( F a \{f}U{u},F < ,k,bg ), 

S" := (FaU{b},F <t k,g). 

* Otherwise, if there is no such h, then: 

• Redehne S := (S \ {5}) U {5", Si ,..., Sj} where with the 
notation of Construction |2.6| (v): 

S' := (F^,F A: + l,^), 

Sj := (^F^,F^,k+ l,ga^j . 

— Otherwise, if there is no such /, then: 

* Redehne S := (S\ {S'}) U {S'} where S' := (Fa,F<j,k + 1, g). 

• Dehne S' := S. 

Output: S'. Then S' is a triangle mush that is equivalent to S and consists of 
triangular systems. 


Proof. Note that we use only operations described in Construction |2.6| for in¬ 
stance, the replacement of S by (S \ {S}) U {S', S"} is an application of, hrst, 
operation (i) and then operation (ii). Therefore, S' is equivalent to S. As each 
S := (Fa,F K \,k,g) £ S' satishes k = r, each element of Fa is constant, i.e., S is 
triangular. 

It remains to show that Algorithm |2.7| terminates. To this end, consider the 
infinite digraph G' = (V ', E') where V' is the set of all semi-triangular systems 
over K[T 1 ; ...,T r ] and, given vertices S\,S 2 £ V, the edge (Si,S 2 ) £ E' exists if 
and only if Algorithm |2.7| replaces Si within a single iteration of the while-loop by 
a triangle mush S' with S 2 £ S'. Let G = (V,E) be the subgraph induced by all 
semi-triangular systems that are reachable by a path starting in S. 

Consider a path (Si,S 2 , ...) in G, i.e., Si £ V and (Si, Si + 1 ) £ E for all i. We 
write Si = (J^, F'^ki, gf). By construction, ki < ki + 1 < r holds for all i. This 
means there is i\ £ Z>i such that fcj+i = ki for all i> i\ and Algorithm 2.7 will 


perform the polynomial division V f = ah+u , i.e., operation (ii) of Construction 2.6 
for each such Si. Since always deg Tfc (6) = 0 holds, we have deg Tfe (/) 
deg Tfc + i(w) and the reduction step only reduces Xb + i-degrees, the sequence 


> 


(Ni)^, Ni := ^ deg Tk ^ 


,(/) e 


^>0 
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is monotonically decreasing. As Ni £ Z>o holds, this sequence either is finite or 
becomes stationary. Assume the latter holds, i.e., there is i 2 £ Z >, (1 such that 
Ni = N i+ i is valid for all i > i 2 . This implies, that for all * > i 2 in the polynomial 
division step only the “&-part” will be added, i.e., 

^ = FuU {&}• 

In particular, the ideal (LT(.?□)) is contained in (LT^^ - )) for each i > i 2 . As 
IK[Zj,..., T r \ is noetherian, the chain 

(LT(JS)) C <LT(^ +1 )) C ... 

becomes stationary, i.e., there is *3 £ Z>i such that (LT(Jq)) = (LT^^ - )) holds 
for all i > 13 . Moreover, as b = LCfc+i(/i) holds and h £ red^^), we have 

LT( 6 ) i (LT (red (•?□))) 2 <LT(^)>. 

Then b cannot be an element of for * > * 3 , a contradiction. Thus, the sequence 
(Ni)i is finite. In turn, this forces the (Si, S 2 ,...) to be finite and acyclic. 

Since each vertex S £ V is adjacent to only finitely many vertices, the previous 
argument shows that G is a finite tree. In particular, the while-loop in Algorithm |2.7| 
will be executed at most |G| times for each vertex S £ V, i.e., the algorithm 
terminates. □ 


Remark 2 . 8 . Algorithm |2.7| is similar to the decomposition into simple systems 
used in (3] . Note, however, that they are interested in special properties (e.g., dis¬ 
jointness) of this decomposition whereas ours is weaker but needs not use operations 
like gcd or subresultant computations. 


An example computation with Algorithm |2.7| will be performed at the end of the 


next section in Example 4.3 


3. Solvability 


We now come to steps (ii) and (iii) in the list on page [I] as before, we as¬ 
sume we are given an ideal I = (/i,...,/ s ) C K[Xj,..., T r \ and a polynomial 
g€K[Ti,..., T r \ and want to answer the question whether there is a: € IK satisfy¬ 
ing l[lj). 

Using Algorithm |2.7| of the previous section with input / and g , we obtain an 
equivalent triangle mush <S that consists of triangular systems. Note that we can 
replace each system (J] 3 >F<a,k,g) £ S with = {0} by the equivalent system 
(0, systems with Fu nK* / 0 clearly are not solvable. Then 0 can be 

rephrased as the question, whether there is x £ IK such that 


f(x) = 0 for all f £ F <t g{x) ± 0 

holds for some (0, F<j, k, g) £ S. Consequently, it suffices to present methods for 
the case S = {5} of a single triangular system. Here is an overview of the steps to 
test whether V ( S ) 7 ^ 0 holds: 


K[Ti ,... ,T r \ 

Ul 

Si- 


rrn 


L[T fel ,... ,Tfc s ] 

ui 

-l(S) h 




Lpfcj,..., TfcJ Lpfcj 

Ul 

*- (0i -7<], k', g'\ 1!- 


/(K) 

UJ 


dense 


dense, monic 


min. polyn. monomial? 


Here, L is a suitable function field. The following proposition reduces the treat¬ 
ment of a triangular system in K[Tj,..., T r ] to a triangular, dense system in 
L [Tfcj,... jTfc,], i.e., a triangular system ( 0 , {/ 1 ,..., / s }, k, g) such that the set 
{fci, • • •, k s } coincides with {fc(/i),..., k(f s )}. 
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Proposition 3.1 (Swap the field). Consider a triangular system S := (0, fc,y) 
in K[Ti,..., T r \. Write I r <| = {/i,...,/ s } and let ki := k{fi ) £ Z>i be as in 
Definition \2.1\ Under the canonical embedding 

l: K[7i, ..., T r \ —¥ L[Tfe 1 ,...,TfcJ, L := K(T); i {ki,... ,k s }) 

we obtain a triangular system l(S) := (0, t(.F<), s, i(g)) that is dense in the polyno¬ 
mial ring L [T *^,... , T^J. Moreover, we have 

V{S) ^ 0 «=> V{i{S)) ± 0. 


For the proof of Proposition 3.1 we recall from m Ch. VI] the generalization of 
evaluation homomorphisms; we will need this to control the elements in L. A place 

K with a subring R £ CL such that 


is a IK-homomorphism e : R v 
x £ 1L. \ R £ 


i 


£ R e and e{x ) = 0. 


Given x £ K r s , denote by e' x : K[T, : ; i £ {ki ,..., fc s }] -A IK the evaluation 
homomorphism. According to ns Thm. 5 in VI.4], we have 


K[T i; i<£ {k u ...,k s }] c R ex c L 



with a place e x : R Sx —> IK extending e’ x . Moreover, we define the domain of t = 
(t \,..., t s ) £ L as the intersection 

S 

Dom(t) := P|Dom(ti), Dom(tj) := jy € K r s ; t % £ R £y j . 
i=1 

Lemma 3.2. In the situation of Proposition \3. 1\ assume we have k\ = 1,..., k s = 
s. Then the following claims hold. 

(i) Consider x £ K and t\,...,t n £ L satisfying E x (t\ ■ ■ ■t n ) = 0. Then 
there is 1 < j < n such that £ x (tj) = 0 . 

(ii) For each t £ V(i(S)) C L* and each x £ Dorn(f) C ]K r s , we have 
( £ xifi)i ■ ■ ■ > e x(t s )i x ) € V(S) where the closure is taken in K r . 

(iii) Given x £ V(S) C K , write x = (x",x’) with x' £ IK , x" £ IK . Then 
there is t £ V{l{S)) C L such that 

x' £ Dom(t) C K r s and (s x >(ti),... ,s x >(t s )) = x”. 

Proof. For (i), we relabel t\,...,t n such that there is k £ Z > 0 with U £ R Ex for all 
i < k and t t (f R Sx for i > k. By definition of places, E x {t~ l ) = 0 for al \ i > k and 
thus 

k / n n \ / n \ / n \ 

n^) = n ^) = e * ( ) ( n ^ & 1 )) = °- 

i—1 \i=l i=k -\-1 / \i=l / \i=k +1 / 

For (ii), given / £ \J (J r <) : g, we have t(/) £ ^ (t(J r <,)) : i(g), which means 
i(f)(t) = 0. Write / = £„ a v TC From 

s r 

f(£ x (t 1 ),...,£ x (t s ),x) = n x 7-s = £ xW)(t)) = 0 

v j=l j=s+l 

we infer that (£ x (t \),..., £ x (t e ), x) £ IK r is an element of the closure V(S) = 
■ 9 ) in K r . 
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We come to (iii). We first show by (finite) induction on 0 < to < s, that there 
are t m + 1 ,..., t s £ L such that for the evaluation homomorphism 

6m '■ K[Ti,..., T r ] —* L[Ti,..., T m ], Tj i —> 

we have (/ m +i,..., / s ) C ker(0 m ) and £ x '(tj) = Xj holds for each m < j < s. 
Nothing is to prove for m = s. Assume now that this claim holds for a fixed 
1 < m < s; we show that it also holds for m — 1. Since we have LC m (/,;) | g, 
g{x) 0 and £ x '(tj ) = Xj for m < j < s, setting a := LC m (/ m ), we obtain 

£x' (0 m (a)) = a(x) | £ x (g) = g(x) ^ 0. 

In particular, 9 m (a) ^ 0. Therefore, the non-zero univariate polynomial f' m := 
6m(fm) £ L[T m ] can be decomposed into linear factors 

n 

f'm = - t mj ) with tmj £ L, c £ L*. 

3 =1 

Note that c = 9 m (a) holds and thus £ x '(c) ^ 0. Moreover, using again £ x '{tj) = Xj 
for j > m and f' m = 6 m (f m ), we have £ x > (f' m {x m )) = f m {x) = 0 where the vanishing 
is due to a; £ V(S). The identity 

/ n 

0 — — Ex' I E (.Xm t 7n j ) 

V i=i 

together with statement (i) provide us with 1 < j < n such that £ x '(t m j) = 
x m . Defining t m := t m j , the elements tm.,... ,t s £ L satisfy the claims: we have 
(fmi ■ ■ •; fs ) — ker (9m— i) since 9 m — i(/m) = = ^ and £x'(trn) — x m holds. 

Using this argument, we now have a map 9q such that both (T<j) C ker((?o) and 
£x'(tm) = x m hold. Setting t := (ti,... ,t s ), we obtain 

t £ U(.F 4 )\U( ff ) = V(l(S)) C L 

because f m (t) = 9 0 (t m ) = 0 for each 1 < m < s and e x ’(Oo(t)) = g(x) ^ 0 implies 
in particular that 9 0 (t) = g(t) ^ 0. By construction, £ x ’(t) = x" holds. □ 


tj, m < j < s, 
Tj , else 



holds, clearly Dom (tj) = K r s \ U(l) is non-empty. If tj ^ 0, we consider the 
product / of the minimal polynomial of tj over L with its common denominator 
and thereby obtain a polynomial h: 

m m 

f = ]T«W\ h = £ K [Tj-, j £ {k u ..., k s }\ [X] 

i=0 i—0 


where h(tj x ) = tj m f(tj) = 0. By definition, each x £ K r s with x ^ Dom(fj) 
must satisfy E^tJ 1 ) = 0. For all i > 0, from a m -i £ KfT,-; j ^ {ki,... ,k s }] we 
know that a m -i £ R Sjs holds and therefore obtain £ x (a m _itJ l ) = 0. We have 


Ex(^m) — f h(tj ) ^ ' Qm—itj J — 0, 


i=l 
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from which we infer that a m (x) = 0 and therefore x £ V(a m ) C IK hold; note 
that the inclusion V(a m ) C IK is proper since a m 7 ^ 0. In other words, 

Dom(tj) D IK r s \V(a m ) 7 ^ 0. 

As finite intersection of supersets of non-empty open subsets, also the set Dom(f) = 
Dom(ti) n...fl Dom(f s ) is non-empty; this completes the proof. □ 


For the remainder of this section, we write L for a field as in Proposition |3.1[ 
note, however, that the following claims also hold for any field L. 

The next step is to make all coefficients of a dense triangular system monic. We 
will call a triangular system (0,.F<| ,k,g) in L[Ti,..., T r \ monic if LC k(f){f) = 1 
for all / £ P < . For instance, the system in Example |2.3| is monic. 

Proposition 3.3 (Make monic). Consider a triangular system S := (0, J-^, fc,g) 
in the ring L[Tj,..., T r \ that is dense in L[T„,..., T r \ for a 1 < n < r. Assume 
there is f £ J-<j with k(f) = n such that T := \ {/} is monic. Then the class 

h £ R . := L[T n+ i,... ,T r \/(P) of h := LCt„(/) is annihilated by a polynomial 

p = bX 3 + fX 3+1 £ L[X] \ {0} with f £ L[A], b £ L* 

where j £ Z>o is maximal with X 3 \ p. Moreover, writing f = hT™ + c with 
m £ Z>o and c £ L [T n ,... ,T r \ such that deg Tn (c) < m, we have a monic dense 
triangular system S' that is equivalent to S: 

S' := (<D,TU{f'},k,g) with f :=T n m -®c £ L [T n ,...,T r ]. 

Lemma 3.4. (i) Consider a triangular system S := (0, ^< 3 , A:, < 7 ) in the ring 

Lpi,..., T r \ that is dense in L[T n ,...,T r ] for a 1 < n < r. Setting 
R := L[T„,..., T r ]/ (.F<|), the ring extension L C R is integral. 

(ii) Let L C R be a ring extension, I C R an ideal and h £ R such that 
h £ R/1 is integral over L. Define J := \fl : h C R and let 

p = bX j + fXi +1 £ L[A'] with f £ L[X], b £ L* 

be the minimal polynomial of h where j £ Z>o is maximal with X 3 \ p. 
Then h! := —f(h)/b £ R yields hh! — 1 £ J. 


Proof. For (i), we write = {/ ra ,...,/ r } and assume k(fi) = i. Define Rj := 
L[T n ,..., T r ]/(fj ,..., f r ) for n < j < r and R r +i = L. The canonical projection 

tt: R j+ i[Tj\ -»• Rj = Rj+i[Tj}/(Jj), f / + (/") 

gives us 7r (fj(Tj)) = fj(Tj) = 0. Since Rj = Rj + i[Tj] and fj £ is monic, 

the generator Tj is integral over Rj+i and non-zero. This shows that in the chain 
R = R n £> ... D R r+ i = L each ring extension is integral, and so is R D L. 

We come to (ii). Note that p(h) £ I and I C J ensures p(h + J) = 0 + J. We 
have 

p(h) = ( hf{h)+b)h J = 0 £ R/J. 

Observe that h is not a zero-divisor: for each x £ R with xh £ \fl : h, already 
x £ \fl : h holds. That is hf(h) + b = 0. Setting h' := —f(h)/b, we obtain 
h'h — 1 £ J from 




= 0 £ R/J. 


□ 


Proof of Proposition\3.3[ Note that the system (0, T, k , g) in L[Ti,..., T r \ is dense 


in L[T ra+ i,... ,T r ). By Lemma 3.4 (i), the residue class h £ R is integral over L, 
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i.e., p exists. Using the inclusion of the ideal yj (J 7 ) : h C L[T n+ i,... ,T r \ in the 
ideal y/ (J 7 ) : g C L[Ti,..., T r \, we obtain 

hti - 1 € : h C yf{F) : 9 with ti := £ L[T n+1 ,..., T r ] 

from the second statement of Lemma [374} One directly verifies the equality of ideals 

■ 9 = V(?) + (/) : 9 = V(F) + </') ■■ 9- 

In particular, V(S) = V(S') holds with the dense triangular system S'. Moreover, 
LCi(/ / ) = 1 by choice of f and S' is monic. □ 

In order to make Proposition |3.3| computational, we first show how one can 
compute the required minimal polynomials. 

Algorithm 3.5 (MinimalPolynomial). Input: an element g £ R where L C R is 
an integral ring extension of finite dimension d := dirndl?). 

• Choosing a suitable L-vector space basis of R, we consider M := [g°,... ,g d ) 
as a d x (d + 1) matrix over L. 

• Compute the kernel K := ker(M) 7 ^ {0}. 

• Choose q £ K C L d+1 such that max(l < j < d; qj 7 ^ 0) is minimal. 

• Define p g := q 0 X° + ... + q d X d £ L[X]. 

Output: p g £ L[A']. This is the minimal polynomial of g £ R. 

Proof. By construction, we have p{g) = Mq = 0. For the minimality, let p' = 
9'j^ e L[X] be the minimal polynomial of g. Then Mq' = Yfj=o9'j^ = 
p'{h) = 0, i.e., q' £ K. By choice of q , we have 

deg(p') = max(l < j < d\ q'j 7 ^ 0 ) > max(l < j < d; qj 7 ^ 0 ) = deg(p). □ 

Remark 3.6. In Algorithm |3. 5 [ the element q £ K can be computed using Gauss¬ 
ian elimination. 


Algorithm 3.7 (Make monic). Input: a triangular system S := (0, k, g) that 
is dense in L[Xi,... ,T r \. We assume = {/ 1 ,..., f r } with fc(/;) = i. 

• For n = r down to 1, do: 

— Set := {/); i > n} C L[T n+ i ,... ,T r ] and define the dense trian¬ 
gular system ( 0 , J 7 ^, k, g). 

- Decompose f n = hT d + c with d £ Z>i and h £ L[T n+1 ,... ,T r ], 
c £ L[T„,..., T r \ such that deg T ^ (c) < d. 

— Use Algorithm |3.5| to compute the monic minimal polynomial ph £ 
L[A] of7)eL[T„ + 1 ,...,T r ]/(J-™). 

— Decompose ph = bXi +1 + aX J with b £ L[A], a £ L* by choosing 
j £ Z>o maximal with X° \ p^. 

— Define b! := — b(h)/a. This yields hb! — 1 £ yj(Xff) : h. 

- Redefine /„ as T d +h'c £ L[T n ,..., T r \. Then S' := (0, ^U{/„}, k, g) 
is a monic triangular system that is dense in L[T n ,..., T r \. 

Output: S'. Then S' is a monic triangular system that is dense in L[Ti,... ,T r ] 
and is equivalent to S. 


Proof. Note that the minim al p olynomial ph exists by Lemma 3.4 (i) since the 
system is dense. By Lemma 3.4 (ii), h £ L[T„ + i,... ,T r \/y/(Fff) : h. is invertible. 
The remaining steps are correct by Proposition |3.3| □ 


We now show that the existence of solutions of a monic, dense triangular system 
can be tested by determining a minimal polynomial. 
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Proposition 3.8 (Solvability). Let S := (0,J r <| ,k,g) be a monic triangular system 
that is dense in L[Ti,... ,T r \. Set R := L[Ti,..., T r ]/(J r <|). Then L C R is an 
integral extension and with the minimal polynomial p g G L[A'] of the residue class 
g G R we have 

V(S) ^ 0 <*=> p g G L[A'] is not a monomial. 

Lemma 3.9. In the situation of Proposition^^ letp G L[A] be a polynomial with 
p(g) £ vW- Then there is k G Z>o such that p g \p k . 

Proof. By assumption, there is k G Z>i such that p(g) k G (.F<|), i.e., p k (g) = 0 G R. 
The monic greatest common denominator a := gcd {p k ,p g ) G L[A] satishes /(<?) = 
0 G R since p k (fg) = p g (]I) = 0. By minimality of p g , we obtain p g = a \ p k . □ 

Proof of Proposition \3.8[ Given x G C L r , the corresponding evaluation 

homomorphism e x fits into the commutative diagram 



The fact, that L C R is integral is Lemma 3.4 (i). Assume now p g = X n holds 


for some n G Z>o, i.e., g G R is nilpotent. By the diagram, g(x) = p x {g) then also 
is nilpotent for each x G V(J : < 1 ) C L . This means g(x) = 0. 

For the reverse direction, assume g{x) = 0 holds for each x G C L , i.e., 

by the diagram, we have p'{g{x))) = 0 with p' := X G L[A]. By Lemma 3.9 there 
is k G Z> 0 such that p g \ (p') k = X k . □ 


We now put the previous propositions and algorithms together to obtain an 
algorithm to check the existence of solutions of a triangular system. This completes 
steps (ii) and (iii) of the list on page [l] 

Algorithm 3.10 (IsSolvable). Input: a triangular system S = (J-b, J- < , k, g) in 
the ring K[Ti,..., T r ], 

• If Pa n K* is non-empty, then: 

— return false. 

• Consider the triangular system l(S) that is dense in L[Tfc i; ... ,T ]~J as in 
Proposition |3.1[ 

• Use Algorithm |3 . 7| with input l(S) to obtain a monic, dense and equivalent 

system S' = (0, k', g’) in L[T fel ,... ,T ks \. 

• Use Algorith m |3.5| to determine the minimal polynomial p g i G L[A] of the 
residue class g' G L [T kl ,... ,T k J/ (^^i). 

• If p g i is a monomial, then: 

— return false. 

• return true. 

Output: true if V ( S ) ^ 0 and false otherwise. 

Proof. By Proposition |3.1[ A lgorithm |3.5| and Algorithm |3.7| S' is equivalent, monic 
and dense. Proposition |3.8| delivers the stated solvability criterion. □ 


4. Monomial containment test and efficiency 

Putting together steps (i)-(iii) listed on pagejl] we are now able to test whether a 
given ideal I C K[Ti,..., T r \ contains some monomial T v , v G Z> 0 . Afterwards, we 
explore the experimental running time of the second author’s implementation [ 12 ] of 
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the algorithm in perl on a series of random polynomials and compare it with Buch- 
berger’s algorithm. Moreover, we compare its efficiency on the examples polsys50 
from m to algorithms listed in [3J Tab. 1]. 


Algorithm 4.1 (ContainsMonomial). Input: generators f\,.... f s for an ideal 
I C K[Ti ,... ,T r \. 






Define the semi-triangular system S := (J r n,0,O,g) where g 
and Jh := {/i,...,/ s }. 


:= T, • • • T r 


Let S be the output of Algorithm 2.7 applied to {S}. 
For each S £ S, do: 

— If Algorithm |3.10| returns true , then 
* Return false. 

Return true. 


Output: true if £ I for some g £ Z> 0 . Returns false otherwise. 

Remark 4.2. In the second line of Algorithm |4.1| it is more efficient to modify 
Algorithm |2.7| such that it checks for solutions immediately after determining a 
new semi-triangular system. 


Example 4.3. In the setting of Example |2.5| 
mark 4.2 to test whether the ideal I := (/i 
mial. 


we apply Algorithm |4.1| with Re- 
f 2 ) C K[Ti,..., T 4 ] contains a mono- 
we apply Algorithm |2.7| to the triangle mush Sq. It will first 


To this end 

choose the polynomial division for (/, h) := (/ 1 , Zb) to obtain 

T4.f1 = (T 2 -T 3 )T 2 f 2 -u, u = (T| - T 3 T 2 2 )T 4 . 


This yields a new triangle mush Si := {S', S''} where S' := ({/ 2 , u}, 0,0, gTi) and 
S" ({/1 , Z *2 ; T 4 }, 0,0, gf In the next step, we obtain triangle mushes 


5 2 := {({«}, {/ 2 },l,r 4 0 ), ({/ 1 , / 2 , T 4 }, 0, 0, g)} , 

5 3 := {(0,{/ 2 ,«},4,T 45 ), ({/!,/ 2 ,T 4 }, 0,O, 5 )}. 


Algorithm 3.10 verifies that the zero-set V (/ 2 , u)\ V{Tig) is empty by the follow¬ 
ing steps: first, Algorithm 3.7 with input (0, {/ 2 , it}, 4, T±g) will return the monic 
system 


(0,{/2,/ 3 },4,T 4 2), f 3 := (T 2 -T 3 )T 2 2 . 

As k(f 2 ) = 1 and fc(/ 3 ) = 2, we set L := K(T 3 ,T 4 ) and the ring R := L[Ti,T 2 ]/(/ 2 ,/ 3 ) 

_ _ 2 

is integral over L with L-basis (1,T 2 ,T 2 ). We have 

%5 = (T 3 - T 2 )T 2 T 3 T 4 £ R, J\g 2 = (T 3 - T 2 ) 2 T 2 T 2 T 2 = 0 £ R, 


By Proposition |3.8| the algorithm may remove this triangular set, i.e., it remains 
to consider 


S 4 := {({/i,/ 2 ,T 4 },0,O, 5 )}. 

The reduction step will remove the redundant equation f 2 . The next steps provides 
us with 


S 5 := {(0,{/i,T 4 },4, W , 2), ({/ 1 ,T 4 , W , },0,O,2)} 1 u := (T 2 - T 3 )T 2 . 


By Algorithm 3.10 the system S := (0, {/ 1 , T 4 }, 4, u'g) has a solution: similar to 
before, Algorithm |3.7| returns the monic system 

(0,{/ 4 ,T 4 },4X2), U := Tl -T 3 


with k{fi) = 1 and k{Ti) = 4. Setting L := K(T 2 ,T 3 ), the ring extension L C R := 
L[T 4 ,T 4 ]/(/ 4 ,T 4 ) is integral with L-basis (1). Since 

fTg = (T 2 - T 3 )T 1 T 2 T 3 = (T 2 - T 3 )T 2 T| £ R 
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is non-zero, its minimal polynomial p = X — (Ti — T^TfTf £ L[X] is not a 
monomial, i.e., V(S) ^ 0 by Proposition 3.8 Thus, So has a solution as we already 
witnessed in Example |2.5[ In particular, I contains no monomial, i.e., the algorithm 
returns false. 


The remainder of this note is devoted to experimental running times. We ap¬ 
ply the perl implementation |12j of Algorithm 4.1 to a series of random ideals 
(/i,..., f s ) C K[Ti,..., T r \ for fixed 2 < s < 5 and running 1 < r < 10. Moreover, 
setting T \= {/i,..., /,}, we distinguish the cases V(T) = 0 and V{T) ^ 0. 

To make the experimental running times better comparable to Buchberger’s 
Grobner basis algorithm j6j, we have reimplemented the latter in perl in two vari¬ 
ants: the first one is the classical version whereas the second one stops as soon 
as a monomial could be found. Both algorithms as well as the testing sets T are 
available at m- The following graphics show the averages over the successful tests. 




V(T) A 0, 3 = 2 V(T) + 0, s = 3 





V(T) = 0, s = 2 V{T) = 0, s = 3 




V{T) = 0, s = 4 V(T) = 0, s = 5 

On the given set of polynomials, Algorithm |4.1| seems to be competitive when 
V(T) ^ 0 whereas, for V(T) = 0, the classical Buchberger’s algorithm usually 
needs less time. 
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Additionally, we have applied Algorithm |4.1| to the set of examples polsys50 
from ESI; its running time as well as the number of performed additions on a 2.66 
GHz machine with time bound 300 seconds and at most 1 GB of RAM is listed in 
the left-hand side part of the following table. We write “n/a” if the computation 
was unsuccessful either due to time reasons or because it was out of memory. 

Moreover, in the right-hand part of the table, we list some of the running times 
listed in [3, Table 1] on the same examples. We want to stress the fact that the two 
sides of this table are only marginally comparable: not only is the goal different 
m deduces more information on the solutions whereas we test the existence of 
solutions), also the machines and maximal running times / memory are different. 


no. 

time|4.1| 

result 

add.s 

time RC1 

time DW1 

time ATI 

1 

> 300 

n/a 

n/a 

3.5 

0.4 

3.0 

2 

> 300 

n/a 

n/a 

7.4 

7.6 

7.1 

3 

30.89 

1 

4956 

> 3h 

985.7 

7538.0 

4 

> 300 

n/a 

n/a 

> 4 GB 

> 4 GB 

0.2 

5 

0.62 

0 

2449 


6 

2.25 

1 

4239 

0.4 

0.1 

0.2 

7 

> 1 GB 

n/a 

n/a 

> 3h 

7352.6 

> 4 GB 

8 

0.14 

1 

214 


9 

11.75 

1 

10149 


10 

0.21 

1 

517 


11 

0.17 

1 

361 


12 

0.74 

1 

1909 

0.5 

0.3 

0.4 

13 

0.15 

1 

214 


14 

0.23 

1 

442 

0.5 

> 3h 

1.5 

15 

29.82 

1 

6655 


16 

> 300 

n/a 

n/a 

0.9 

1.4 

1.8 

17 

> 300 

n/a 

n/a 

6.5 

4.7 

75.5 

18 

2.34 

1 

4324 

0.3 

0.1 

0.1 

19 

> 300 

n/a 

n/a 

419.9 

0.4 

0.4 

20 

0.25 

1 

668 


21 

> 300 

n/a 

n/a 

1.6 

86.6 

4.5 

22 

> 300 

n/a 

n/a 

0.6 

1.2 

1.5 

23 

> 300 

n/a 

n/a 

0.4 

0.1 

29.5 

24 

> 300 

n/a 

n/a 

1.2 

1.3 

1.0 

25 

0.25 

1 

537 

1.2 

> 3h 

> 4 GB 

26 

1.07 

0 

4610 


27 

1.47 

1 

2320 


28 

9.89 

1 

6632 


29 

0.15 

1 

297 

0,3 

0,3 

0,3 

30 

> 300 

n/a 

n/a 

> 4 GB 

> 4 GB 

45.3 

31 

> 1 GB 

n/a 

n/a 

> 4 GB 

> 4 GB 

> 3h 

32 

0.41 

1 

1200 


33 

> 1 GB 

n/a 

n/a 

3.4 

1.3 

3.5 

34 

> 300 

n/a 

n/a 

911.5 

> 3h 

> 4 GB 

35 

> 300 

n/a 

n/a 

1.5 

1.2 

1.7 

36 

0.13 

1 

160 


37 

0.27 

1 

633 


38 

11.10 

1 

6878 


39 

> 300 

n/a 

n/a 

0.6 

1.2 

0.6 

40 

> 300 

n/a 

n/a 


41 

> 300 

n/a 

n/a 

1.5 

1.5 

7.0 

42 

0.35 

1 

1028 


43 

> 1 GB 

n/a 

n/a 

0.7 

3.1 

0.2 

44 

> 300 

n/a 

n/a 

24.5 

3.4 

1.2 

45 

> 300 

n/a 

n/a 


46 

> 300 

n/a 

n/a 


47 

16.40 

1 

10465 

1.3 

2.8 

13.0 

48 

0.23 

1 

563 


49 

> 300 

n/a 

n/a 

0.3 

610.2 

0.5 
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